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An unfitted /^-interface penalty finite element 
method for elliptic interface problems 

Haijun Wu* Yuanming Xiao^ 

Abstract 

An hp version of interface penalty finite element method (hp-IPFEM) 
is proposed for elliptic interface problems in two and three dimensions on 
unfitted meshes. Error estimates in broken H 1 norm, which are optimal 
with respect to h and suboptimal with respect to p by half an order of p, 
are derived. Both symmetric and non-symmetric IPFEM are considered. 
Error estimates in L 2 norm are proved by the duality argument. 
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1 Introduction 

£> ' 

_~ ■ Let £! = £7iUrufi2bea bounded and convex polygonal or polyhedral domain 

^^ | in M. d , d = 2 or 3, where J7i and O2 are two subdomains of $1 and T = dtti C)d£l2 

OO ■ is a C 2 -smooth interface (see Fig. [1} . Consider the following elliptic interface 

Cn ' problem: 

Q ■ ( - V • (a(x)Vu) = /, in fii U D, 2 , 

(1) \[u}=9d, [(a(a;)Vu) • n] = g N , on T, 

u = 0, on dfl, 

k> , where [v] — w|o! —v\n 2 denotes the jump of v across the interface T, n is the unit 

5^ ■ outward normal to the boundary of Oi, and a(x) is bounded from below and 

above by some positive constants. Note that a{x) is allowed to be discontinuous 
across the interface T. 

It is well known that the finite difference method (FDM) compared to the 
finite element method (FEM) is easy to implement since it discretizes the PDEs 
on simple meshes. But it is not easy to construct high order FDMs when the 
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Figure 1: A sample domain 51 and an unfitted mesh. 



geometries are complicated. By contrast, the FEM provides a systematic way 
to design high order schemes for problems with complicated geometries, while 
the classical FEMs require that the interfaces and boundaries are fitted well 
enough by underlying meshes in order to achieve optimal convergence rates (cf. 
[121 EI] ) ■ Therefore it has always been attractive to develop systematic ways 
for constructing high order numerical methods on simple unfitted meshes. Un- 
fitted mesh techniques are particularly convenient for time-dependent problems 
when moving interfaces are involved. One could just use the same mesh on the 
domain for different times instead of repeatedly remeshing the domain to fit the 
moving interfaces. Another motivation for developing unfitted methods is from 
the adaptive finite element methods. We know that in an adaptive finite ele- 
ment procedure, the meshes are locally refined repeatedly to equidistribute the 
error. But it is readily to encounter the "inverted elements", i.e., the elements 
with negative directional areas/volumes (cf. |37j). near curved interfaces during 
interface-fitted local refinements, in particular, in the three dimensional case. 

Various finite difference schemes have been proposed on unfitted Cartesian 
grids, for example, the immersed boundary method by Peskin [34], the ghost 
fluid method by Osher and his coworkers [T7], the immersed interface method by 
LeVeque and Li [26 , and the matched interface and boundary method by Zhou, 
et al. 41 . These methods have the first order, second order or higher order 
truncation errors near the interfaces. They have been generalized in various 
ways and applied to many problems [30l EH [38l E21 [24l EH [29l SOI etc]. 

The existing finite element schemes on unfitted meshes are usually designed 
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by proper modifications of the standard finite element methods near the inter- 
faces. The immersed finite element method by Li, et al. 28, 19 modifies the 
basis functions for nodal points near the interfaces to satisfy the homogeneous 
jump conditions. The linear immersed finite element method is optimally con- 
vergent in H 1 and L 2 norms for two dimensional problems [13] . We refer to 
|llj for adaptive immersed interface finite element methods for solving elliptic 
and Maxwell interface problems with singularities. The multiscale finite element 
method by Chu, Graham, and Hou |14) modifies the basis functions near the 
interface by solving appropriately designed subgrid problems. The method was 
proved, for the two dimensional case, to be of first-order accuracy in H 1 norm 
and of second-order accuracy in L 2 norm. Note that when the interface inter- 
sects an element in a straight line, the multiscale finite element method coincides 
with the linear immersed finite element method (cf. |14j). The penalty finite 
element method proposed by Babuska [4J modifies the bilinear form near the 
interface by penalizing the jump of the solution value across the interface. The 
method may use high order elements but is suboptimally convergent in h. The 
linear case of the method was analyzed again by Barrett and Elliott [7] and was 
proved to be optimally convergent in h under more regularity assumptions on 
the exact solution. We remark that the idea of adding penalty terms is widely 
used in the interior penalty Galerkin methods [21 El [TBI HH1 E3 etc] . The un- 
fitted finite element method proposed by A. Hansbo and P. Hansbo [2 3) can be 
viewed as an improvement of the linear version of the Babuska's method. This 
unfitted method was proved, for two dimensional elliptic interface problems, to 
be optimally /i-convergent in both H 1 and L 2 norms and the error estimates 
are uniform with respect to the relative position of the interface to the unfitted 
grid. The key idea of their method is using weighted average flux across the 
interface in the bilinear form instead of using the arithmetic average one as in 
the standard interior penalty discontinuous Galerkin methods. In their method, 
the weights along the intersection of the interface and a grid element are cho- 
sen according to the areas of the two parts of the element separated by the 
interface. More recently, Massjung [33 proposed an hp-unfitted discontinuous 
Galerkin method for Problem (p} which uses weighted average flux across the 
interface and penalizes not only the jump condition on the solution value but 
also the jump condition on the flux. It was proved that, for the two dimen- 
sional case, the method converges in broken H 1 norm at an optimal rate with 
respect to h and at a suboptimal rate with respect to p by a factor of p. The 
flux penalty term is not essential in the convergence analysis but helpful for 
numerical stability (cf. [33]). We also refer to [51 [35] for recent more studies on 
the discontinuous Galerkin methods for elliptic interface problems. 

In this paper we propose an /ip-interface penalty finite element method (hp- 
IPFEM) for the interface problem ([I]). In this method we penalize the jump 
conditions on the solution value as well as the flux across the interface and adopt 
the arithmetic average flux across the interface in the bilinear form. We provide 
a rigorous analysis to show that, for both two and three dimensional cases, the 
error estimates in broken H 1 norm are optimal with respect to h and suboptimal 
with respect to p by a half order of p. Both symmetric and non-symmetric 
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IPFEM are considered. Some /ip-error estimates in L 2 norm are also derived 
by the duality argument. We would like to mention that our tricks on dealing 
with discontinuities across the interface are common in the standard interior 
penalty discontinuous Galerkin methods. They are adopted by [331 G2I] except 
the arithmetic average one. Sure we do not add penalty terms on edges/faces of 
elements away from the interface as the standard interior penalty discontinuous 
Galerkin methods do. 

For any intersection segment/patch e between the interface and a grid el- 
ement K e , let Kf — K e n fli,i = 1,2, be the two parts of K e separated by 
the interface. The key idea of our analysis is the use of the following straight 
forward identity on e: 

— = Vi H — (U1-U2), where Vi = (dVv h )\K? • n, i = l,2. 

By this we may bound the term of the arithmetic average flux (i.e. 1 ' 1 +'" 2 ) by 
any component in the average (say w,-, i = either 1 or 2) and the jump of flux 
(i.e. v\ — V2) which is controlled by the flux penalty term. This component is 
then estimated by a local inverse trace inequality on the corresponding part of 
the interface element. As a result, it suffices to prove the local inverse trace 
inequality on either Kf or iff instead of on both of them as treated in [33] . 
Therefore we avoid discussing various relative positions of the interface to an 
element K e which are already complicated in two dimensional case since Kf 
or K% may not be shape regular. This makes our analysis much simpler and 
easy to apply to the three dimensional case. To the best of our knowledge, our 
results give the first hp a priori error estimates for interface problems in three 
dimensions and the best in two dimensions. 

The rest of our paper is organized as follows. We formulate the /ip-interface 
penalty finite element methods in Section [2] and list some preliminary lem- 
mas in Section |3j The H 1 - and L 2 - error estimates of the symmetric interface 
penalty finite element methods are given in Section|4l The error estimates of the 
non-symmetric interface penalty finite element methods are given in Section O 
Section [5] is devoted to the proofs of the local trace inequality and the local 
inverse trace inequality. 

Throughout the paper, C is used to denote a generic positive constant which 
is independent of h, p, and the penalty parameters. We also use the shorthand 
notation A < B and B > A for the inequality A < CB and B > CA. A ~ B is 
for the statement A < B and B < A. Note that the constants in this paper may 
depend on the jump of the coefficient a(x) across the interface. Although it is 
an interesting topic to derive error estimates with explicit dependence on the 
jump of the coefficient a(x), it may bury our basic idea in a more complicated 
analysis. We leave this issue to future work. 
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2 Formulation of /^-interface penalty finite ele- 
ment methods 

To formulate our /ip-IPFEMs, we need to introduce some notation. The space, 
norm and inner product notation used in this paper are all standard, we refer 
to [9j [15] for their precise definitions. 

Let {Mh} be a family of conforming, quasi- uniform, and regular parti- 
tions of the domain il (cf. [9]) into triangles and parallelograms/tetrahedrons 
and parallelepipeds. For any K G Mh, we define hx '■= dia,m(K). Let 
h := maxKeMh hx- Then Hk ~ h. Note that any element K G Mh is consid- 
ered as closed. Obviously, each partition Mh induces a partition, denoted by 
Ih , of the interface T (cf. Fig. [I| . Note that any interface segment/patch e G Xh 
is either contained entirely in one element in Mh and has nonempty intersec- 
tion with its interior, or is the common edge/face of two neighboring elements 
in Mh- For any e elj, let K e G Mh be one of the element(s) containing e and 
let Kf = K e f\Vti,i = 1,2. 

Introduce the "energy" space 

(2) 
V := {v : v\ Ut = v t , where v t G flft(fi), v t \ K G H 2 (K), MK G M h , i = 1,2} . 

For any K G Mh, let V P (K) denote the set of all polynomials whose degrees in 
all variables (total degrees) do not exceed p if K is a triangle/tetrahedron, and 
the set of all polynomials whose degrees in each variable (separate degrees) < p 
if K is a parallelogram/parallelepiped. Denote by U% the /ip-continuous finite 
element space, that is, 

(3) Ul := {v h 6 #o(^) : »*M P P W, VK G M fc } . 

Wc define our interface penalty finite element approximation space VF as 

(4) Vjf := {v h : u/,1^ = u»ft, where v ih G E/£, i = 1, 2} . 

Clearly, V? G V G L 2 (fi). But V^ ^ II 1 (0). 

We also define the jump [v] and average {v} of u on the interface T as 

(5) [w] :=u|ni-u|n a , {>} := ^ • 



Testing the elliptic problem (fT]) by any dgF, using integration by parts, and 
using the identity [vw] = {v} [w] + [v] {w}, we obtain 

\] / aV« • Vw — / {aV« • n} [v] — gN {v} = / fv 
i=1 Jo; Jr Jr Jn 
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(6) ct,h(u,v) :=^2 aVu ■ Vv - J^ / ( {aVu • n} [v] + /? [u] {aVu ■ n} 

»=1 "' - e(El h "' e 

+ J (u,v) + Ji{u,v), 

(7) J„(v):=El^/H[i 

e£l h K ° Je 

(8) J x (u, v):=J2 11 ^r 1 [ I aVu • n l [« Vw ' n l > 

where j3 is a real number, 70 and 71 are nonnegative numbers to be specified 
later. Define further the linear form Fh(-) on V: 

(9) F h (v) := [ fv+ [ g N {v} - p f g D {a\7v • n} + J D (v) + J N {v), 

Jn Jr Jt 

„2 



(10) J D ( v ):=J2lT- [9d[v], 

(11) jjf^^^TJ^L [ gN [aVvn] 



It is easy to check that the solution u to the problem (p} satisfies the following 
formulation: 

(12) a h {u,v) = F h (v), \fveV. 

We are now ready to define our /ip-interface penalty finite element methods 
inspired by the formulation (fT2"|) : Find Uh £ V£ such that 

(13) a h {u hl v h ) = F h {v h ), Mv h e VJf. 

Remark 2.1. (i) The above tricks on dealing with the discontinuities are 
from the interior penalty discontinuous or continuous Galerkin methods 
(see, e.g., JH [751 \TU\j). When j3 — 1, ah(-,-) is symmetric, that is, 
dh(v,w) — ah(w,v), the method is termed as the symmetric interface 
penalty finite element method (SIPFEM), which corresponds to the sym- 
metric interior penalty Galerkin method [2, 391. On the other hand, 
when j3 ^ 1, ah{-,-) is non-symmetric. In particular, when f3 = —1. 
the method is called the non-symmetric interface penalty finite element 
method (NIPFEM), which would correspond to the non-symmetric inte- 
rior penalty Galerkin method studied in fSSjj. In this paper, for the ease of 
presentation, we only consider the case (3 = 1 and fi = — 1, nevertheless 
the ideas of this paper also apply to the case (3 =/= ±1. 

(ii) The unfitted finite element method given in j23f does not include the 
penalty term J\ and uses a weighted average {v} — Kiv\q 1 + K2v\n 2 in- 
stead of the arithmetic one in (|5|), where Ki\ e = \Kf\ / \K\. Optimal error 
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estimates in both H l and L 2 norms were derived in [231 for the linear 
element (p — 1) in two dimensions. 

(iii) Massjung [33] proposed an hp-unfitted discontinuous Galerkin method for 
Problem ([lj which uses weighted averages similar as in J23f and penalizes 
both the jump of the solution value and that of the flux. It was proved that, 
for the two dimensional case, the method converges in broken H 1 norm at 
an optimal rate with respect to h and at a suboptimal rate with respect to 
p by a factor of p. 

Next, we introduce the following broken (semi-)norms on the space V: 



£*(«<) *-f h K " 



£^ll[ a V,.n]||! Y 

„CT. " / 



(15) IIMIU := ( \Hl h + E ^ IK aVu ' n >l' 2 



3 Some preliminary lemmas 

Recall that, for any interface segment/patch e E Ih, K e £ M-h is an element 
containing e and Kf — K e fl fli, i = 1,2. We have the following lemma, which 
gives the trace and inverse trace inequalities on K e , whose proof will be given 
in Section El 

Lemma 3.1. There exists a positive constant ho depending only on the interface 
r and the shape regularity of the meshes, such that for all h £ (0, ho] and any 
interface segment/patch e E Ih- the following estimates hold for either i e = 1 
or i e = 2: 

(i) For anyv E H l {Kf e ), 

(16) IMI i2(e) < h K \' 2 \\v\\ L2{KV + \\v\\%\ KV \Wv\\ l ' 2 \ KV . 

(ii) For any v h E V p (K e ), 

(I 7 ) \\Vh\\ L 2 (e) ^7172 \\Vh\\ L 2 {K! )• 

We remark that, if e G dK e , then the result (i) in the above lemma is a 
direct consequence of the standard local trace inequality, and the estimate in 
(ii) is already known (cf. [10]). 

The error analysis relies on the following well-known hp approximation prop- 
erties (cf. HUG!!]): 
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Lemma 3.2. Suppose w £ H s (O) fl Hg(fl),s > 2. Let jj, = min{p + 1, s}. 

(i) There exists a piecewise polynomial Wh, wh\k G T^p(K), such that, for 
j = 0, I,--- ,s, 

ft* 4- -? 

Ik - «>/.||ff:f(K) ^ "^TJ Ikllff.(Jf) , ^ K g Mi- 



(ii) There exists Wh G ?/£ smc/i iftai 

_ffere i/ie invisible constants in the above two inequalities depend on s but are 
independent of h and p. U^ is the hp-continuous finite element space defined in 
©• 

The following lemma is from the standard Sobolev extension property (cf. 

my 

Lemma 3.3. Let s > 2. There exist two extension operators E\ : H s (p,\) \— > 

H s (n)nH^(n) an dE 2 ■. {w <= H s (a 2 ) ■■ w\ dn = 0} ^ H s (n)nH^(n) such that 

(E i w)\n i =w and ||-E»Hlff«(n) i$ MLs^n*) > * = 1 . 2 , 
where w € H s (Qi) for i = \ and w € H s (fl 2 ), w\on = for i — 2 (cf. Fig. QJJ. 

4 Error estimates for the symmetric interface 
penalty finite element method 

In this section we derive the H 1 - and L 2 - error estimates for the /ip-SIPFEM 
(i.e. the case ft = 1). The following lemma gives the continuity and coercivity 
of the bilinear form a h (, •) for the SIPFEM. 

Lemma 4.1. We have 

(18) M«,itf)|<2||H|| lfh ||H|| lihl Vv,w£V. 

For any < 71 < 1 and < s < 1, there exists a constant ao, £ independent of h, 
p. and the penalty parameters such that, if Jo > ao,e/jii then, for < h < ho, 

(19) a h (v h ,v h )>(l-e)\\\v h \\\l h , V^e^, 
where the constant ho is from Lemma \3.1i 
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Proof. (fT5|) is a direct consequence of the definitions ([S])-©, (fT4")) -(fll) j) , and the 
Cauchy-Schwarz inequality. It remains to prove (|19[) . We have, 



2 

(20) a h (v h , v h ) = V U 1/2 V^ - 2 V / {oV% • n} [ Vfc ] 

i=l y '' e£l h Je 



7o P 



l\h K ' 



p 



^E 

Mli.ft - E — ^IKaV^-n 



V/iJHLa(e) + "^2" lli aV ^ • n ]llz,2( e , 

hx? „ r _ -, ,,2 



eei fl 



7o V 



l£ a (e) 



- 2 53 / {aVt; ft • n} [w h ] . 



It is clear that, for any eGl;,, 



(21) 



{aVwft ■ nf | e = (oVvfc)!^-? ■ n ■ 



(-1)' 



[aVw/j • n] | e , 



where z e is specified in Lemma 13. II From Lemma l3.1l (ii). there exists a constant 
C\ such that, 



(22) (aVth)|Kj -n 



< 



P 



P 






O^Vtffc 



£ 2 (tffJ 



Therefore, from flU) and ([22]), 

( 23 ) E^H{ aV ^' n }ll^ 



eGXh 



7o P" 



(e) 



/(/ 



^ 2 E ^5( ||(« v ^)i^ • n L^ + i ii[« v ^ • n ]ii^ { . 



eel, 



7o P 



1 

\L 2 (e) ' 4 



< V — 

t^ 7o 

e£lh 



2C? a 1 ' 2 ^ 2 



L*{Kf 



E^ii[« v ^- n ]i^ 



eGlh 



270 P 2 



(e) 



r 2C ? x ^ 11 ii=2 

<maxl ,- J «/, Nft- 
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On the other hand, 

(24) 2^ f {aVv h ■ n} [v h ] 
eei h Je 

- 2 E (f( aV OK ' n 

a 1/2 Vv h 



eeih 

< E (^ 

eex h v 



+ 2 llt aVw ' 1 ' n ]lli 2 (e) J II Wlli»(e) 



P 



vvcnW" L JllL2(e) 



< 



E 



lot 



aVuft -n]|| L2(e) ||K]|| i2(e 



4Cl 2 i/2r, 2 



L *( K i) 4 fee 



II Ml 



2 



1 /lA'e I, r _ ,||2 . £ "foP" |i i|,2 



70 £ p* 

MC? 1 e\ . .a 
< max I , , - % j h . 

V7o£ 7o7i^ 2 ' 



K' 



By combining flgOJ), ([23]). and (|24[). we conclude that ([HI) holds. This completes 
the proof of the lemma. D 

The following lemma is an analogue of the Cea's lemma for finite clement 
methods. 

Lemma 4.2. There exists a constant ao independent of h, p, and the penalty 
parameters such that for < 71 < 1, 70 > ao/71, and < h < ho, the following 
error estimate holds. 

\\\ u - u '»llli h ~ inf ll u ~~ ^llli h ' 

where the constant h is from Lemma \3.1l 

Proof. From (fTSj) and (fT5]) . we have the following Galerkin orthogonality 

(25) ah(u-Uh,Vh) = 0, W/, e V/f. 

For any z^ E V^ 1 , let 7^ = u^ — Zh- By letting e = 1/2 and in Lemma [4.11 
denoting by ao = a ,i/2, an d using (|2"5)) . we have, for < 71 < 1 and 70 > 
ao/71 > 

Hhhllli,/, <2ah(r)h,T)h) = 2a h (u - Zh,f]h) < 4 |||u - «h||| lih HWHi,/, 1 

and hence, 

\\\uh - z h \\\^ h <4:\\\u-Zh\\\ lth , 
which implies that 

III" _ U h\\\l,h = III" ~ Z h + z h - Uh\\\i th < 5 |||U - Z h \\\ lh . 

This completes the proof of the lemma. □ 
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The following lemma gives some approximation properties of the space VF. 

Lemma 4.3. Let s > 2 be an integer and let (i = min{p+ 1, s}. Suppose the 
solution to the elliptic interface problem (Q]) satisfies ufa € H s (fli),i = 1,2. 
Then there exists iih £ V£ such that, for < h < ho, 



(26) \\u-u h \\ L2{n) <—(J2\\ u \\%(n t 



P -. 



/l 1 7l\!/2 h^ 1 /v^ 9 x ■'■ ' 

(27) u - & h \ lih < (- + 70 + — + ^) 4=372 E U B.«l, 

\p 70P p J p s ■ i l A \^— ■? y * 



where the constant ho is from Lemma \3.1[ 

Proof. Let it.; = -EiU, where _Ei,i = 1,2 are the extension operators from 
Lemma 13.31 Then 



(28) Ui GH s (Q)nH^n), Uik=«k ) Ikll^n^Ni^)' 
From Lemma T3. 2f ii) and (|28|) . there exist Uj/j e E/F, i = 1,2 such that 

h^-i 

(29) Ik - Uih\\ H i(ci) ~ ^_J ll u llff*(n 4 ) > J = 0,1. 

Define w ft <E V£ by M/Js^ = w^lr^i = 1,2. It is clear that (j2"Bf holds. 

Denote by ( = m - % and by £j = itj — Uih,i — 1,2. Obviously, £|n; : 
Ci|ni)« = 1,2. Next we estimate each term in |||C||li /,. (cf. (|14| - (IT5)) '). Clearly, 

<*» EIK'HIU,, 5 ?^^? 11 " 11 "^)- 

From Lemma 1331 (i) and (12U1) . 



(=i 






£ E F^E fe iic<ii V.) + n^n^(ice) iivcj l2( ^ 

S ^E ( ft_1 IK*lli»(n) + IKJl 2 (o) IIVCIIl^o) 

To estimate the remaining terms, we estimate ||(aV£)|n; ■ n llz,2( e ) , i = 1,2, e S 
Zh. From Lemma l3.2l (i). there exist piecewise polynomials &»/,, Ujfc|if G V P (K), 
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such that 

h^~ j 
\\ui -v.ih\\ H i(K) ~ ~^j \\ u ^\\h-(k) i Vj = 0,1,2, i = 1,2, if eMft. 

Therefore from Lemma |3. 11 

||(aVC)k • n||^ 2(e) < ||V(ui - Wih)||i2 (e) + ||V(w ifc - u ih )f L2{e) 
<h K l \\V(ui -Uih) f L 2 (Ke) 

o 
P 2 

+ T ||V(Ujfe - Uj + Uj - U»h) || r,2fK-e1 

2 , -P 



< 



h2fj,-3 ,;l 



n 2s-4 



If 



* II #«(«■«) + ~T \\V( U i ~ U ifc)llz,2(#e) 



It follows from ([28) and O that 



7,2^-3 

ElKaVOk-nlli-w^^zrllull^,. 



Thus, 



( 32 ) E f^ll[ aV C-n]|li 2(e) + ^ll{«VC-n}||i 2(e) 



2 



4 1+ ^(j>M 



Then (J2TJ) follows by combining ([30 ]) -(|32 ]) , This completes the proof of the 
lemma. D 

By combining Lemma 14.21 and Lemma 14.31 we have the following theorem 
which gives error estimate in the norm 1 1 1 - 1 1 1 x h for the SIPFEM. The proof is 
straightforward and is omitted. 

Theorem 4.4. Let s > 2 be an integer and let /i = min {p + 1, s}. Suppose the 
solution to the elliptic interface problem ([T]) satisfies w|q. G H s (Qi),i = 1,2. 
Then there exists a constant ag independent of h, p, and the penalty parameters 
such that for < 71 < 1 and 70 > cto/'fi, the following error estimate holds: 

lll»-«fclHi,h^7o 3=372(53 H tt llH«(n«) ) » V0</i<h o , 



where the constant ho from Lemma \3.1\ depends only on the interface Y and the 
shape regularity of the meshes. 
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Next we derive the L 2 -error estimate by using the Nitsche's duality argument 
(cf. [15]). Consider the following auxiliary problem: 

{— V ■ (b(i)Vw) = u — Uh, in fii U 0,2, 

[w] = 0, [(o(i)Vto) • n] = 0, on T, 

w = 0, on <9fi, 

It can be shown that w satisfies (cf. [4]) 

(34) IMIh2 (0i) + |MI H 2 ( n 2) < ||« - u h\\ L 2 {n) ■ 

Theorem 4.5. Under the conditions of Theorem \4-4\ there holds the following 
estimate for the SIPFEM. 

W ( 2 2 V /2 

where the constant h is from Lemma \3.1l 

Proof. Let r\ = u — U] x . Testing (J53")) by r\ and using (f^5|) we get 

(35) hllL^O) = a h( w i V) = a h(v, u>) = a h(v, w~w h ), 
where Wh £ V£ satisfies the estimate (cf. Lemma |4~3|) 

(36) \\\w - w h \\\ lh < ll' 2 -rjl Y^ IMIff2 (a) . 

1 i=\ 



Therefore from (fl81) and (134 

-, „||| < ^,1/2 

1/ 



v\\ 2 L 2 { n) < 2 Ilk " «fc|lli,h IIMIIi,/, £ to Zlfi IMIi^n) IIMIIi,/, - 



that is ||?y|| i 2/ ) < 7 ^72 Uhllli /,, which completes the proof of Theorem 14.5 
by using Theorem 14.41 □ 



5 Error estimates for the non-symmetric inter- 
face penalty finite element methods 

In this section we derive the H 1 - and L 2 - error estimates for the /ip-NIPFEM 
(i.e. the case (5 = —1). The following lemma gives the continuity and coercivity 
of the bilinear form a h (-, •) for the NIPFEM. 

Lemma 5.1. We have, for the NIPFEM ((3 = -1), 

(37) M«.«')l<2||H|| ljh ||Hlli,h. Vv,w£V. 

(38) ah(vh,Vh)> 7 TV IIMIIih) Vvh£VF, V0<h<h , 

7o7i + max (71,1) 

where the constant ho is from Lemma \3.1l 
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Proof. (|3"T)) is a direct consequence of the definitions (JS])-©, ([T3|) - (fT5|) . and the 
Cauchy-Schwarz inequality. From (|23[) . 



Ilklll?,, < (l + max (_i,__))|| Wh ||; 

v v 7o ^7o7i 7/ 

Then (|38[) follows from the identity ah(vh,Vh) — ||f/i||i «,- This completes the 
proof of the lemma. □ 

Following the argument of Theorem 14.41 we have the following error estimate 
in the norm 1 1 1 - Ilia, h ^ or tne ^P-NIPFEM. The proof is omitted. 

Theorem 5.2. Let s > 2 be an integer and let fj, = min{p + 1, s}. Suppose the 
solution to the elliptic interface problem ([T]) satisfies u^ G H s (tti),i = 1,2. 
For < 7i < 1 and jq > l/7i, the following error estimate holds. 



^ - 1/2 



where the constant h is from Lemma \3.1l 

The following theorem gives the L 2 -error estimate for the /ip-NIPFEM. 

Theorem 5.3. Under the conditions of Theorem \5.2\ there holds the following 
estimate for the NIPFEM. 

. fofJ, ^/j-1/2 / 2 \ 1/2 

h- u h\\i^ (m ^{lfo^ + -r^)\^h\\H'[a i )J > V0<h<h Q , 

where the constant h is from Lemma \3.1l 

Proof. The proof is similar to that of Theorem 14.51 So we only sketch the 
differences here. (1551 becomes 



(39) IMI^n) = a h(w, rj) =ah{v, W )~ 2 Y1 / M ^ aVw ' n ^ 

eex„ Je 

=a h (ri, w - w h ) - 2 ^ M (aVw ■ n) , 
eei„ Je 

where r\ = u — Uh, W is defined in (|33[) , and u>h £ V? satisfies the estimate in 
HMD. From dSZD,®, and ®, 

(40) a h (ri, w~w h )< ll' 2 ~rj^ IMIz^fi) \\\v\\\i,h > 
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The last term in (I39[) is bounded by using the Cauchy-Schwarz inequality and 
the trace inequality: 

(41) - 2 E / W ( aVw • n ) ^ 2 E II W H^(e) lloVw • n|| L2(e) 

^^^(e^uvhiW)' 172 

<-^J-Jo(»?,»7) 1/2 ||Vt£;|| ia(r) 

7o P 

<^-J ( VlV )^(\\Ww\\ L2{ni) \\Vw\\ H1{ni V L> 

7o P 

h 1 / 2 

^-T^-INIIi,hll''lli a (n)' 
7o P 



where we have used (|34[) to derive the last inequality. By combining (|39|) - (j4T 
we conclude that 

/ 1/2 h h 1 / 2 \ 

\\v\\ma)Zho -rr 2 + -772- ) IIMIU . 

v P ' 7 ' p y 



which completes the proof of Theorem 15.31 by using Theorem 15.21 □ 

Remark 5.1. Here is a brief comparison between the SIPFEM and the NIPFEM: 

• The stiffness matrix of the SIPFEM is symmetric while that of the NIPFEM 
is not. 

• Both methods require that the penalty parameters satisfy 7071 > C but 
C is a problem dependent constant for the SIPFEM while C can be an 
arbitrary positive constant for the NIPFEM. 

• The L 2 -error estimate for the SIPFEM is optimal with respect to h while 
that for the NIPFEM is not. 

6 Proof of Lemma 13.11 

In this section we prove the local trace inequality and the local inverse trace 
inequality in Lemma 13.11 The proof is divided into two cases, the two dimen- 
sional case and the three dimensional case. Recall that, e € 1/, is an interface 
segment/patch, K e € J\4h is an element containing e, and Kf — K e DH,i, i = 1,2. 

6.1 Two dimensional case 

Since T is a C 2 boundary of fii, it can be expressed as a union of open sub- 
segments Ti,i £ A := {1, 2, • • • , N} such that each I\ is parameterized by some 
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function x = r<(£) e [C 2 (ii)] 2 satisfying |r^(£)| ^ 0, where Ii is an open inter- 
val. Since I\ n Vj is either empty or an open subsegment of T, e is contained 
entirely in some I\. if /ix= is small enough, and hence e is parameterized by 
x = r(£),£ 6 / with r(£) = rj.(£) and / C Ik- It is clear that the length of I 

<h K e. 



strip of width Pk' 




P 

Figure 2: An interface segment e and the element K e in 2D 

Let Po be any fixed point on e and let Tp be the tangent line to T at Pq. 
Let P be a point in if 6 which achieves the maximum distance, denoted by h, 
from the tangent line Tp . Denote by pk c the diameter of the largest disk 
contained in K e . Clearly, we have h > Pk<= /2 > /i_k-<= (see Fig. [2] otherwise i^ e 
would be completely contained in a strip of width less than px c along the axis 
Tp , which causes a contradiction). Now, we may choose Kf to be the one of 
two subregions Kf and iff which contains the point P (in Fig. [2j for instance, 
if 4 e e = Kf). Without loss of generality, we assume P is the origin. Suppose Po 

corresponds to the parameter £ = £o- Then the vector PPo = r(£o)- It is clear 
that 

(42) \r(£o)\<h K ., |r%)|~l. 
We have the following Taylor's expansions at £o, 

(43) r(0 = r(£ ) + r%)(£ - Co) + 0(h 2 K .), 

(44) r'(0=r'(^ )+O(^e). 

To proceed, we need to bound the following term G(£) from below: 

|r(£)xr'(0| 



G(0 



|r'(OI 



where a x b = ai&2 — 02^1 denotes the scalar cross product of two dimensional 

vectors. Recall that |axb| = |a||b| sin (a, b). Note that G(£ ) = * '/ \ 

l r (Co) 
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equals to the distance from P to Tp , thus 

Gfa) = h>h K .. 
By combining (|4"2")) - (|4"4")) . we derive that 



Hto)xr'(Zo)+Q(h? K .)\ 
\r'(Z )+O(h K s)\ 



Q(0= ' ,L/^ ■ ^M. \ ' = G(fr) + 0(fee). 



Then there exists a positive constant /io depending only on the interface F and 
the shape regularity of K e , such that if < fee < /i , then 

(45) G($>h K *, for all f e I. 

Next we prove the trace inequality ([TBI . Let if 6 be a fan-like region in Kf , 
with eU PD U PE a,s its boundary (cf. Fig. El, that is, 

K e := {x : x = i • r(£), t G [0, 1], £ G /}. 

Let v G C l (K£ ) and consider its restriction on e as follows, 

(46) V 2 (r(0) = J ^((t 2 v 2 (t-v(0))dt=f 2{t 2 vv t +tv 2 )dt. 
From vt = Vv • r, we have 

v 2 (r(0) < I 2{\v\\Vv\\r\+v 2 )tdt. 
Jo 

Integrating along e, we find 
(47)Jv 2 ds = y>(r(£))|r'(£)|d£ 

< 2 f f {\v\\\7v\\r\+v 2 )t\r'(£_)\dtd^ 

= 2jJ (| w ||V«||r|+ V 2 )-t|r(Oxr'(C)|-^ydtde 

where we have used |r(£)| < fee to derive the last inequality. By noting that 
£|r(£) x r'(£)| is the absolute value of the Jacobian determinant of the mapping 
x = t ■ r(£) and using (|45l, we obtain that 

(48) |M| 2 L2(e) < . f * / (\v\\Vv\h K e +v 2 )dx 

ml ?e /G(£) fee 

< IH^K^Vvll^^ + h K l\\v\\ 2 L2(liar 
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Since K e C Kf , we conclude that 

— 1 1 1 

IMU 2 (e) < h K l\\v\\ L 2 (K e e ) + \\v\\l HK ej\\Vv\\l 2{K! j. 

By a density argument, the above inequality is also valid for v 6 H (K? ). This 
completes the proof of p^|) . 

It remains to prove the inverse trace inequality (|17[) . We need the following 
inverse inequality on V p ([0, 1]) (cf. [55]): 

(49) IKIb([o,i])<P 2 |KIU2([o,i]) Vto h e^p([0,l]). 

By noting that Vh(t ■ r(£)) is a polynomial in £ of degree < p if if e is a triangle 
and of degree < 2p if i"f e is a parallelogram, we have from the above inverse 
inequality that 

(50) 

vl(r(0)= [ 1 ^(t\l(t-r(0))dt = 2 [\v h (t-v(0)^(tv h (t-r(0))dt 



<2^\tv h (t-r(0))y (£ (§- t (tv h (t-r(0))) 2 



<p 2 / ^ 2 (i-r(£))di. 



o 
Following the same guidelines of (j47]l and (j48|). we derive that 

(51) \\v h \\h {e) = J^Mmr^Ol^Zp 2 Jj\l(t-r(0)t\r'(t)\dtdt; 

„2 



< P 



f P 

/ v h dx ^T \\ v h\\ 2 L 2 ( Kt 



which yields the conclusion of (|1T|) . 

6.2 Three dimensional case 

The argument for the two dimensional case is readily extended to the three 
dimensional one. We only sketch the proof by indicating the necessary modifi- 
cations in the proof for the two dimensional case. 

Since T is a C 2 interface, it can be expressed as a union of open subpatches 
I\, i e A :— {1, 2, • • • , N} such that each Ti is parameterized by some function 
X = rj(£, rj) G [C 2 (Ui)] 3 satisfying r^ x r iri \ ^ 0, where U t is an open domain 
in K 2 . Since I\ (~l Tj is either empty or an open subpatch of T, e is contained 
entirely in some I\. if hx" is small enough, and hence e is parameterized by 
X = r(S,,ij), (£,r]) E U with r(£, rt) — !"&(£, 77) and U C Uk- It is clear that the 
diameter of U < h,K<> ■ 

Let Pq be any fixed point on e and let Tp be the tangent plane to T at Pq. 
Then there exists a point P € K e such that the distance from P to Tp > hj(* ■ 
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:tangcnt plane 
to Tk at Pq 



Figure 3: An interface patch e and the element K e in 3D 



Without loss of generality, we assume P is the origin. Suppose Pq corresponds 
to the parameter (£,77) = (£o,??o)- We have 

(52) |r(£ ,»7o)| < h K ; M&, Vo) X r,,(£o,?7o)| ~ 1. 

(53) rfer?) - r(£ 0) »/o) + r^ , %)(£ - Co) + r„(£ ,»to)fa - »») + 0(J&.)> 

(54) r £ (f,?y) x ^(£,77) = r e (£ ,r7o) x r»,(£o,?7o) + 0(h K *)- 



Let 



G&77) 



■fo^-M^ 7 ?) xr^(^,77))| 



Mf,»?) xr,((,r;)| 

By noting that G(£o:?7o) is the distance from P to Tp , we have G(£o,r)o) > 
fix?. From (l52j) - (|54|) . there exists a positive constant fto depending only on the 
interface T and the shape regularity of K e , such that if < hie < ft.O) then 



(55) 



G(t,v)Zh K . V(^,»?)eJ7. 



Next we prove the trace inequality (fTB")). Let X e be the cone with apex P 
and curved base e as shown in Fig.[3J i.e., 



K e :={x: x=t-r(S,n), iG[0,l], K,»?)ei7}. 



We have 



1 a 



(56) v 2 (r(H,v)) = J Q- t (t 3 v 2 (t-r(t;,n))dt 



(2t 6 vv t + 3t 2 v 2 ) dt 



Therefore, from |r| < hx' and the fact that t 2 |r(^, 77) • (^(^,77) x v v (£,,v))\ is 
the absolute value of the Jacobian determinant of the mapping x = t ■ r(£,7y), 
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we conclude that 

(57) 

v 2 do = / w 2 (r(£,7?))|r ? (£,77) Xr,((,i])| d£d?? 



< 



f f (\v\\Vv\\T\+tP)f > \T^1 1 )XT v (^,Tl)\dt^dv 

Ju Jo 



<j J (| 1 ,||VH^+ l , 2 )t 2 |r(e,r ? ).(r 5 (C,77)xr,(e,r ? ))|^-ydidedry 
<—f 1 rut \ l {\v\\Vv\h K e+v 2 )dx 

<lkllL 2 (^)l|V«|| L2(if e ) +^|M| 2 2(i?c) 



which implies that (|T6|) holds. 

It remains to prove (fl"7|) . By following the proof of ([50)) , we have 



vl(r(Z,v))<P 2 I t 2 v 2 (t-r^,rj))dt. 
Jo 

Then, from (55J), 

(58) /i£d<7 = / t£(r(£,r,)) |r e (e,»|) X r„(&»,)| d^dr, 

<p 2 [ [ v 2 (t-r^,rj))t 2 iTtfoTfixTnfoffildtiitdT, 

Ju Jo 

o 



^)L^ dS; ~fc IKI1 ^)' 



which implies that (IT71) holds. This completes the proof of Lemma \'3. II 
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